Infrared properties of propagators in Landau-gauge pure Yang-Mills theory at finite 

temperature 
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The finite-temperature beiiavior of gluon and of Faddeev-Popov-giiost propagators is investigated 
for pure SU{2) Yang-Mills theory in Landau gauge. We present nonperturbative results, obtained 
using lattice simulations and Dyson-Schwinger equations. Possible limitations of these two ap- 
proaches, such as finite-volume effects and truncation artifacts, are extensively discussed. Both 
methods suggest a very different temperature dependence for the magnetic sector when compared 
to the electric one. In particular, a clear thermodynamic transition seems to affect only the electric 
sector. These results imply in particular the confinement of transverse gluons at all temperatures 
and they can be understood inside the framework of the so-called Gribov-Zwanziger scenario of 
confinement. 

PACS numbers: ll.lO.Wx 12.38.Aw 12.38.Gc 12.38.Lg 12.38.Mh 14.70.Dj 25.75.Nq 



I. INTRODUCTION 

Thermodynamic observables, such as the free energy, 
show a thermodynamic transition in Yang-Mills theory 
at a finite temperature Tc P, This phase transi- 
tion separates a low-temperature phase — which is ex- 
pected to be highly non-perturbative and characterized 
by quark and gluon confinement — from a (in princi- 
ple) perturbative high-temperature phase, where color 
charges should be deconfined. Indeed, due to asymp- 
totic freedom, the running coupling constant g{T) van- 
ishes with increasing temperature [3] and the quark-gluon 
plasma can be viewed as a weakly-coupled system. At the 
same time, across the phase transition, the long-distance 
potential between fundamental test-charges changes from 
a linearly-rising confining potential to an essentially flat 
one 01, i.e. the (temporal) string tension is zero at high 
temperature. 

On the other hand, non-perturbative phenomena have 
to be expected even for an arbitrarily small coupling at 
the renormalization scale In particular, several stud- 
ies have already pointed out that non-perturbative ef- 
fects should be present in Yang-Mills theory at any tem- 
perature, i.e. the high-temperature phase should also be 
highly non-trivial. For example, it is well-known that the 
magnetic sector is affected by strong non-pcrturbative in- 
frared (IR) problems [IjQ- Moreover, at finite temper- 
ature, the spatial string tension is non- vanishing Q. Fi- 
nally, the formal infinite-temperature limit of the theory, 
i.e. the 3-dimensional reduced theory 0, @, is still a con- 
fining theory 3, [M [H [H, [Tsj . Thus, a clear description 
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of Yang-Mills theory as a function of the temperature T is 
still lacking and the fate of confinement at large T is still 
an open question. In particular, one should reconcile the 
confining properties of the dimensionally-reduced theory 
with the vanishing of the conventional string tension at 
large temperature. 

In order to make contact with the continuum formu- 
lation of SU{Nc) Yang-Mills theory, it is necessary to 
consider quantities living either in the su{Nc) alg ebra or 
in the continuum gauge group SU{Nc)/Z]s[ [Ij]. This 
clearly limits the usefulness of center observables such as 
the Polyakov line, often employed in lattice simulations 
and in effective theories. On the other hand, propagators 
of the elementary degrees of freedom of the theory, such 
as gluon and ghost fields, fulfill the above condition and 
provide access to non-perturbative aspects of the theory 

MM- 

In Landau gauge, at zero temperature, the framework 
of the Gribov-Zwanziger [ll|, [13] and of the Kugo-Ojima 
[m scenarios provides a basis for understanding the con- 
finement mechanism in Yang-Mills theory. In particular, 
these scenarios predict that the ghost propagator should 
be IR enhanced, compared to the propagator of a mass- 
less particle. At the same time, the gluon propagator 
should vanish in the IR limit. These predictions are 
suppor ted by results obtained u sing different methods 

Mm m m m m m m, m, mm m mmm- 

In addition, several calculations indicate that the pre- 
dictions are also (at least partially) valid in the high- 
temperature phase ItI. Isil. [3^ [ssl and in the infinite- 
temperature limit Bniiniiiii,!!, m, [M m, m Hi] . 

Here we evaluate gluon and ghost propagators for the 
SU(2) case using Landau gauge, extending earlier numer- 
ical studies that focused on the high-temperature phase 
[m and on the infinite-temperature limit U [13, S, HB] ■ 
We present nonperturbative results, obtained using lat- 
tice simulations and Dyson-Schwinger equations (DSEs). 
In Section |ll] we present a short review of the Lorentz 
structure for the Landau-gauge gluon and Faddeev- 
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Popov-ghost propagators at finite temperature. Results 
using lattice gauge theory for these two propagators at 
various temperatures, both below and above the thermo- 
dynamic transition, will be reported in Section IIIII In 
that section we also discuss the influence of finite- volume 
effects on the numerical results. In Section ITVl the same 
propagators are studied in continuum space-time using 
DSEs, which allows us to consider the limit of vanishing 
momentum p — > 0. Let us recall [Tlj that imposing the 
minimal Landau gauge condition, i.e. restricting to the 
Gribov region, does not modify the DSEs, even though 
it provides supplementary conditions for their solution. 
We believe that the simultaneous use of these two non- 
perturbative methods — lattice simulations and DSEs — 
allows a better understanding of our results and of their 
physical implications. A possible interpretation of the 
results from lattice and from DSEs is presented in Sec- 
tion [V] The main components of this interpretation are 
a spatial/magnetic sector, which is essentially unaffected 
by temperature and by the thermodynamic transition, 
and a temporal/electric sector, which does show a sensi- 
tivity to the transition. Let us stress that the results from 
DSEs suggest that the difference between the two sectors 
occurs for any non-zero temperature and not only above 
the phase transition. In particular, the Gribov- Zwanziger 
and Kugo-Ojima confinement mechanisms remain quali- 
tatively unaltered in the magnetic sector when the tem- 
perature is turned on. Finally, we present our conclusions 
in Section rvTl Some observations about the use of asym- 
metric lattices in numerical simulations are reported in 
Appendix]^ Analytic and numerical details of the DSE 
analysis, presented in Section lTVl are collected in Appen- 
dices E El 

Preliminary results have been presented in [37| . 



II. PROPAGATORS AT FINITE 
TEMPERATURE 

At finite temperature, the Euclidean symmetry is man- 
ifestly broken by the thermal heat bath. As a conse- 
quence, the propagator of the gluon, which is a vector 
particle, can no longer be described by a single tensor 
structure. Instead, one can consider the decomposition 

M 



Dt{p) = PL{p)Dt\p) + PM(p) 



(1) 



where the two tensor structures P^,y{p) and P^,y{p) are 
respectively transverse and longitudinal in the (three- 
dimensional) spatial sub-space. Of course, in the full 4d 
space both tensors must be transverse, due to the Landau 
gauge condition. In terms of the momentum components 
Pf^ we can write (with no summation over repeated in- 
dices) [11] 



P^iP) = (l-^po)(l-'5.o) (<5^.- 

PLip) = Pf^Ap) - pL{p) 



Pl^pl^ 



(2) 
(3) 



Clearly, the tensor P^^^p) is the transverse projector 
in the three-dimensional spatial sub-space, while P^i,{p) 
is the complement of P'^^p) to the four-dimensional 
Euclidean- invariant transverse projector P^v{p) = Sfj,i, — 
PtJ-Pv/p^- At the same time, the decomposition ([T]) de- 
fines the two independent scalar propagators D^{p) and 
Df^{p), which are respectively the coefficients of the 3d- 
transverse and of the 3d-longitudinal projectors. Since 
one expects a different IR behavior for these two func- 
tions, it is not useful to contract the gluon propaga- 
tor with the usual 4d-transverse projector P^i,{p), i.e. 
to consider a linear combination of the two scalar func- 
tions [p] and Df^ (p) . By contracting the propagator 
D'^Kp) with Pifj,{p) or with P^i,{p), inserting a unit ma- 
trix in color space and using the relation 



PoA^o 



(4) 



valid in Landau gauge, we find (for non-zero three- 
momentum p) the relations"'^ 



Drip) 



1 



{d-2)Ng 

3 2 

< J2 K(P)^U-P) - ^ ^oipWoi-p) > (5) 
^1 P 

pI 



D^ip) ^ —\l + ^^< Al{p)Al{-p) > . (6) 

Here, d 4 is the space-time dimension, Ng = 3 is the 
number of gluons [i.e. Ng — — 1 in the SU{Nc) case] 
and the summation over color indices is implied. Clearly, 
for po = the 3d-transverse propagator Dt{p) coincides 
with the 3d gluon propagator, while the 3d-longitudinal 
propagator D^lp) can be identified with the propagator 
of the would-be-Higgs field in the three-dimensionally- 
reduced theory. Thus, we can associate the 3d-transverse 
gluon propagator to the so-called magnetic sector and 
the 3d- longitudinal propagator to the electric sector. Of 
course, at zero temperature the two scalar propagators 
coincide, i.e. Dt{p) = Dl{p) = D{p). 

As for the ghost propagator Dg{p), since it is a scalar 
function, no additional tensor structures arise in this 
case. 

A second consequence of considering a theory at finite 
temperature is that the propagators depend separately 
on the energy po and on the spatial three-momentum \p\ . 
Moreover, the energy is always discrete, i.e. po = 2TiTn 
with n integer. Of course, in studying the IR properties of 
the theory, we are mainly interested in the zero (or soft) 



^ Similar considerations apply also to the gluon propagator when 
studied using asymmetric lattices. In Appendix |X] we present 
results for the gluon tensor structures on a strongly asymmetric 
lattice. 
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TABLE L Data of the configurations considered in our numerical simulations. The three-dimensional case will be discussed 
in Section [III Al The Setup column will be used in the text to refer to the various numerical setups. The value of the lattice 
spacing a has been taken from ^3^] for the three-dimensional case and from jZ^ for the four-dimensional case. Sweeps indicates 
the number of sweeps between two consecutive gauge-fixed measurements. The temperature has been evaluated using Eq. 
([8]). In the column Method we indicate a reference where numerical details (including error determination, etc.) are discussed 
for that particular set of data. With the symbols and * we indicate that (respectively) only the gluon or only the ghost 
propagator has been considered. Note that only the results corresponding to the Setups 16, 17, 19, and 20 have already been 
published. Finally, in the Setup 6 we found one "exceptional" configuration that induced very large statistical fiuctuations on 
the results. Thus, in this case, we needed a much larger statistics in order to acquire the same level of accuracy obtained in 
the other cases. 



Setup 


N*, Nt X Nf or iV^ 


P 


a'^ [GeV] 


Configurations 


Sweeps 


r [MeV] 


K^''^ [fm] 


Method 


1 


12"* 


2.3 


1.193 


426 


40 





1.98 


[36] + 


2 


16* 


2.3 


1.193 


424 


45 





2.64 


[36] + 


3 


20* 


2.3 


1.193 


405 


50 





3.30 


[36] + 


4 


32** 


2.3 


1.193 


30 


320 





5.28 


[24]* 


5 


10 X M"' 


2.3 


1.193 


401 


40 


119 


2.31 


[36] 


u 


1 n V 90'^ 

lU A ZU 




1 1 0"^ 


1 zl9fi 




1 1 Q 


o.ou 


loDJ 


7 


10 X 26^ 


2.3 


1.193 


405 


50 


119 


4.29 


[36] 


8 


4 X 20^ 


2.3 


1.193 


444 


40 


298 


3.30 


[36] 


9 


4 X 26^ 


2.3 


1.193 


405 


45 


298 


4.29 


[36] 


10 


4 X 34^ 


2.3 


1.193 


405 


50 


298 


5.61 


[36] 


11 


3 X 28^ 


2.4 


1.651 


429 


50 


550 


3.34 


[36] 


12 


4 X 40^ 


2.5 


2.309 


360 


50 


577 


3.41 


[36] 


13 


2 X 20^ 


2.3 


1.193 


437 


40 


597 


3.30 


[36] 


14 


2 X 32^ 


2.3 


1.193 


410 


45 


597 


5.28 


[36] 


15 


2 X 42^ 


2.3 


1.193 


400 


50 


597 


6.93 


[36] 


16 


20^ 


4.2 


1.136 


6161 


40 





3.47 


^6] 


17 


30'' 


4.2 


1.136 


10229 


45 





5.20 


[36] 


18 


40'' 


4.2 


1.136 


3306 


50 





6.94 


[36] 


19 


80^ 


4.2 


1.136 


200 


200 





13.9 


[35] + , [41J* 


20 


140'' 


4.2 


1.136 


30 


250 





24.3 


[35]+ 



modes po — 0. The other (hard) modes {po ^ 0) have an 
effective thermal mass of lixTn and seem to behave like 
massive particles 0, [H, [H, |39[ . 

Let us stress that these observations apply both to the 
lattice and to the continuum formulations of Yang-Mills 
theory. In the lattice case one also has to consider the 
definition of the gluon field 

A,(x) = l[[/,(x)-C/t(x)]^_,_ (7) 

as a function of the link variables Ufj_{x), which is based 
on the expansion C/^ (x) 1 + ig^aAf^ (x) . Hence, in the 
high-temperature (symmetry-broken) phase, one should 
consider only configurations for which the trace of the 
Polyakov l oop has a positive value when averaged over 
the lattice [3i|,[43|. 

III. LATTICE RESULTS 

The propagators of the gluon and of the Faddeev- 
Popov ghost have been numerically evaluated using the 
methods described in Refs. [13, Ha, [s^, |4l| . Details can 
also be found in Table [H Finite temperature has been in- 
troduced according to the standard procedure of reducing 



the extent of the time direction compared to the spatial 
ones [i^. We consider Nt lattice sites along the tempo- 
ral direction and Ng sites along the spatial direction with 
Ng ^ Nt. After taking the infinite-spatial- volume limit 
Vs — Ng — > oo , the continuum limit is given hy Nt oo 
and a —^ 0, keeping the product aNt fixed. This yields a 
temperature 

aNt 

In the following, we first consider the symmetric three- 
dimensional case, discussing finite-volume effects in de- 
tail. We then move on to the four-dimensional case and 
discuss the effects of finite temperature. 

A. Finite-volume effects 

Before considering the finite-temperature dependence 
of the propagators, it is worthwhile to discuss the ef- 
fects of finite volume on their IR behavior. Clearly, a 
finite volume V affects all processes with a correlation 
length ^ of the order of (or larger than) V^^"^. [In the 
finite-temperature case one can expect finite-volume ef- 
fects when ^ > Vs^^'''~^\ where Vg is the spatial vol- 
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Gluon propagator ] 




0.5 1 1.5 2 

p [GeV] 

Gluon dressing function | 




FIG. 1: The gluon propagator D{p) (top) and dressing func- 
tion p^D{p) (bottom) as a function of the momentum p for 
various volumes in three dimensions. Open circles, open 
squares, open triangles, full circles and full triangles corre- 
spond to Setups 16 (V = 20^), 17 (V = 30^), 18 {V = 40=^), 
19 {V = 80^) and 20 (V = 140^), respectively. The solid line 
is the result from Dyson-Schwinger calculations jlH, yielding 
a gluon propagator suppressed as p"'^^ in the IR limit. 



ume.] Since confinement is induced by an infinite cor- 
relation length [l5|, finite- volume effects are likely to be 
observed in the numerical evaluation of gluon and ghost 
propagators.^ These effects should be considered when 



comparing the numerical data to the predictions of the 
Gribov-Zwanziger and of the Kugo-Ojima confinement 
scenarios. Let us recall that several calculations in the 
continuum, e.g. using DSEs [111, 

, Il9l | or the renor- 

malization group [22i] , are in favor of these scenarios and 
confirm their predictions. In addition, using DSEs one 
can also study (in the continuum) the consequences of a 
finite volume on these propagators. In particular, it has 
been found [i^ that in a finite volume the gluon propa- 
gator is suppressed at small momenta and finite at zero 
momentum, while the ghost propagator is IR enhanced. 
These results are in qualitativ e ag reement with the find- 
ings using lattice calculations |24| . [27l Isol l46j . Moreover, 
using DSEs it can be argued [43| that an IR-finite gluon 
propagator is indeed an artifact of the finite volume and 
that one obtains a null propagator at p = when the 
infinite-volume limit is considered. On the lattice the 
situation is clearly more complicated. Even in the case 
of a power-de pen dence of the type -D(O) oc , it can be 
very difficult [35| to perform simulations with sufficiently 
large values of the physical lattice volume V so that one 
can really control the extrapolation to infinite volume. 
This is in fact the case and up to now there is not yet 
any convincing result on a (symmetric) four-dimensional 
lattice showing an IR- vanishing gluon propagator D{p). 
The situation is actually even worse than this. Indeed, in 
order to sec an IR-suppressed propagator one should ob- 
tain that D{p) has a maximum for some value of p 7^ 0. 
Up to now this has not been obtained in the Landau 
4d case (using symmetric lattices), even when consid- 
ering very large lattices (see Fig. 10 in [2^ and Fig. 3 
in [30|]). There are, however, evidences of a suppressed 
gluon propagator when considering 4d asymmetric lat- 
tices [4al , even though in this case it is difficult to extract 
quantitative information from the data [13, Hlj. Also, a 
suppressed gluon propagator at small momenta can be 
obtained in the 4d case when simulations are done in the 
strong-coupling regime [4^ . Finally, a suppressed gluon 
propagator has recently been obtained (using symmetric 
lattices and /3 values in the scaling region) for a Landau- 
like gauge condition [i^, confirming that the difficulties 
in finding a similar IR-suppressed gluon propagator in 
Landau gauge are probably related to finite-size effects. 

Numerically, it is clearly easier to consider first the 
three-dimensional case, since one can then use much 
larger lattice sides. Moreover, the Gribov-Zwanziger sce- 
nario can be applied to three dimension s | 49| — as also 
supported by numerical studies [Slj, [H, l35l. [s^ and by 
calculations using DSEs [llj, [l^ — and the IR suppres- 
sion of the gluon propagator is expected to be stronger 
than in the 4d case (i.e. the IR exponent is larger in the 



^ One should recall that finite-size efli'ects are more or less pro- 
nounced depending on the quantity considered. In particular, 
one does not expect such efli'ects for quantities with an intrinsic 
mass-scale m, such as quarks or hadrons, since in these cases the 
relevant correlation length is of the order of l/m. This has an 



important consequence in studies using functional methods: the 
specific behavior in the deep IR re gion of the gluons is nearly 
irrelevant for hadronic observables Il5l . lid |43| and results are 
sensitive only to the behavior at intermediate momenta (about 
0.5-1 GeV). 
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Ghost propagator ] 




0.5 1 1.5 2 2.5 3 3.5 

p [GeV] 



Ghost dressing function | 



2 




FIG. 2: The ghost propagator Dcijp) (top) and the dressing 
function p^Dcip) (bottom) as a function of the momentum 
p for two different volumes in three dimensions. Full trian- 
gles and open circles correspond to Setups 17 {V = 30"^) and 
19 iy — 80^), respectively. The solid line is the result from 
Dyson-Schwinger calculations [l^l, yielding a ghost propaga- 
tor enhanced as p~^'^ in the IR limit. 



3d case). Finally, since the theory is finite, renormaliza- 
tion issues do not obscure the interpretation of the re- 
sults. Let us recall that various lattice calculations have 
been performed in the 3d case, both for the propagators 
[31I [34 [35I HI] and for the 3-point vertices Q . These 
results support the assumptions usually made in DSE 
calculations p^ . 

We now consider symmetric lattices in the 3d case (see 
run parameters at the bottom of Table |T| . These can 
be interpreted as lattices at T = in 3d or as lattices 
at T = cx), simulating the dimensionally-reduced the- 
ory (without the Higgs field). We compare the numer- 



Effective coupling ] 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 



p [GeV] 

FIG. 3: The effective coupling p^ Dl;{p) D{p) as a function of 
the momentum p for Setups 17 (open squares, V = 30^) and 
19 (full triangles, V — 80'^). The solid line is the result from 
the same Dyson-Schwinger calculations [l3| reported in Figs. 
nandU 

ical data to Dyson-Schwinger calculations In Fig. 

[T] we present results in the 3d case for the gluon prop- 
agator using lattice simulations and DSE calculations. 
Considering the lattice data on the larger lattices, there 
is a clear maximum at Pmax ~ 400 MeV. Moreover, as 
the lattice volume increases, the IR suppression becomes 
stronger. However, even using very large volumes (i.e. 
V = 140^ [sEDi only few momenta are available in the 
interval [Q,Pmax\- The comparison to the DSE solution 
(for infinite volume) shows qualitative agreement. Con- 
cerning the gluon dressing function p^D{p), it is found 
that lattice results exhibit a power-law behavior for mo- 
menta below 400 — 500 MeV. The result is very similar 
to the DSE solution, albeit differing by a constant fac- 
tor. Of course, when considering the dressing function 
p^D{p), the factor p^ is dominant in the IR limit and the 
agreement between lattice and DSE results looks better. 
Also note that the pre-factor of the power-law obtained in 
DSE studies is sensitive to the approximations performed 

Results for the ghost propagator Dg[p) and the dress- 
ing function p'^Dg{p) are reported in Fig. [2l In this case 
there are no evident finite-volume effects in the lattice 
data for Dg{p), even though such effects are visible at 
small momenta when considering the dressing function 
P^Dg{p)- In particular, the dressing function shows a 
reduced IR divergence for the largest lattice at small mo- 
menta p. On the other hand, one should recall that, for a 
relatively small statistics, the so-called "exceptional con- 
figurations" are probably not adequately sampled and 
that they contribute significantly to the IR enhancement 
of Dq (p) , i.e. this IR enhancement could be underesti- 
mated [2^, [3^. When considering the comparison to DSE 
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3d-longitudinal gluon propagator at T=0 MeV 



3d-longitudinal gluon propagator at 7=149 MeV 
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3d-transverse gluon propagator at T=0 MeV 



3d-transverse gluon propagator at T=149 MeV 
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FIG. 4: The volume dependence of the gluon propagators in the low-temperature phase. The top panels show the 3d-longitudinal 
propagator while in the bottom panels we report results for the 3d-transverse propagator. The left panels show the results 
at zero temperature {Ns = Nt) and the right panels at T ~ 119 MeV (Nt = 10), well inside the low-temperature phase. We 
always compare results obtained using three different physical volumes. The smallest volumes {Setups 1 and 5) are indicated 
with full circles, the middle-size volumes {Setups 2 and 6) are represented by open triangles and the largest volumes {Setups 3 
and 7) are indicated by full triangles. Momenta are aligned along the x axis. 



results, we see that agreement is again at the qualitative 
level (see in particular the dressing function in the bot- 
tom figure). Indeed, the lattice data show a weaker diver- 
gence than the one obtained using DSE calculations. As 
said above, this could be related to an insufficient statis- 
tics for the ghost propagator when large lattice volumes 
are used. 



Finally, one can consider the effective coupling^ defined 
to be proportional to Dq{p) D{p), for which DSEs pre- 
dict a constant limit at zero momentum. The comparison 
of lattice data and DSE results is made in Fig. [S] It is 



^ Since the 3d theory is (ultraviolet) finite this quantity is of course 
not a running coupling in the sense of the rcnormaUzation group 
as in the 4d case [ISl. |16|| . Nevertheless, using DSE calculations, 
one can make a prediction for its IR, behavior also in the 3d case. 
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3d-longitudinal gluon propagator at T=298 MeV 



3d-longitudinal gluon propagator at T=597 MeV 
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3d-transverse gluon propagator at T=298 MeV 
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FIG. 5: The volume dependence of the gluon propagators in the high-temperature phase. The top panels show the 3d- 
longitudinal propagator while in the bottom panels we report results for the 3d-transverse propagator. The left panels show 
the results at T « 298 MeV {Nt — 4), slightly above the thermodynamic transition, and the right panels at T ~ 597 MeV 
{Nt = 2), well inside the high-temperature phase. We always compare results obtained using three different physical volumes. 
The smallest volumes {Setups 8 and 13) are indicated with full circles, the middle-size volumes {Setups 9 and 14) are represented 
by open triangles and the largest volumes {Setups 10 and 15) are indicated by full triangles. Momenta are aligned along the x 
axis. 



clear that the asymptotic (constant) regime is reached 
in the Dyson-Schwinger resuh only for momenta of the 
order of 200 MeV. For the lattice data, when ^ = 30^, 
we do not have any point in the range [0, 200] MeV and 
it is difficult to say what one could find in the p — 
limit. On the other hand, in the case of the larger lat- 
tice {V = 80^), one sees a coupling decreasing when p 
is smaller than about Pmax ~ 400 MeV. It thus seems 
that the two approaches differ qualitatively in the deep 
IR limit. Nevertheless, it is interesting to note that the 



maximum of the data for the larger lattice is obtained 
for the momentum p = Pmax ~ 400 MeV, where also the 
gluon propagator reaches its maximum value. Thus, a 
possible interpretation [29l . IZsj is that finite-size effects 
are different for the two propagators and one might see 
an IR-finite effective coupling only in the infinite- volume 
limit. In other words, since the prediction from DSE 
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FIG. 6: The /3 dependence of the gluon propagators for a physical spatial volume Vs of approximately (3.35 fm)^ and temperature 
T « 575AdeV, using three different lattice setups. We report data for the 3d-longitudinal propagator in the top-left panel, for 
the 3d-longitudinal dressing function in the top-right panel, for the 3d-transverse propagator in the bottom-left panel and for 
the 3d-transverse dressing function in the bottom-right panel. We represent with full triangles data at P — 2.3 {Setup 13), 
with full circles data obtained ai /3 — 2.4 {Setup 11) and with open triangles results for /3 = 2.5 {Setup 12). Momenta below 
approximately 2.4 GeV are considered along the a;-axis, larger momenta are along the 3d-spatial diagonal. 



studies is based on the relation^ an — '^clg = (4 — rf)/2 
[ill . [23I I. one can imagine that the true values for the 
exponents ao and ac are obtained only in the infinite- 
volume limit and that at finite volume the above relation 
does not need to be satisfied.'^ 



Here we define the IR exponents and using the relations 
D{p) ~ pZ^D-a and Dg{p) ~ p'^'^c-s. 
For a different interpretation, see |50|| . 



Let US note that Dyson-Schwingcr studies find that the 
expected power-law behaviors for the propagators (in 3d 
as well as in 4d) start to appear only for momenta smaller 
than some energy-scale A/ < 200 MeV [i^. Also, one 
should consider with caution lattice results at the small- 
est non-zero momentum 2/asm{Tr/N) w 2T:/{aN), i.e. 
when the quantities considered may "feel" the bound- 
aries. Thus, one should try to extract the IR behavior of 
the propagators only using data in the range 



3in(-^) <p< A/. 



(9) 



9 



With a lattice spacing a of order of 1 GeV~^ one conse- 
quently needs N ^ 50. Present simulations have access 
to this region only in the three-dimensional case. In the 
four-dimensional case the situation is more complicated 
and, in particular, since the IR suppression of the 4d 
gluon propagator is predicted to be weaker than in three 
dimensions (jp-^^ instead of [O, H^, one probably 

needs even larger lattice sides than in the 3d case, mak- 
ing these computational studies very demanding. Look- 
ing at the results obtained in the 3d case and considering 
the data obtained with the largest 4d lattices available 
[29I [30| , one can argue that volumes larger by a factor of 
10 (i.e. a factor of about 2 in lattice side) are probably 
necessary in order to see a decreasing gluon propagator 
also in the 4d case. On the other hand, for the ghost 
propagator and with the available lattice sides, the IR- 
enhancement is clearly observed [2^ . [27l [30l . |40 |. even 
though the value obtained numerically for the IR expo- 
nent og is usually smaller than the DSE result ^11,, 2^ . 



B. Finite temperature 

Here we present our calculations for gluon and ghost 
propagators in 4d at finite temperature. We consider four 
different temperatures, two of them below the thermody- 
namic transition (i.e. at zero temperature and at T « 119 
MeV) and two above it (i.e. one at T « 298 MeV and 
the other at T « 597 MeV). Our runs are summarized in 
Table m Note that, for a time extension Nt = 4, the crit- 
ical coupling is /3 = 2.299 corresponding to Tc « 295 
MeV. Here the calculations for this time extension have 
been done at /3 = 2.3. Thus, the temperature T « 298 
MeV is just above the thermodynamic transition, while 
the highest temperature (T « 597 MeV) corresponds to 
about twice the critical temperature T(.. This allows us 
to make contact with previous studies [3l| of the high- 
temperature phase, which investigated the gluon prop- 
agator for T > 2Tc. In order to clarify possible finite- 
volume effects, we considered for each temperature three 
different spatial volumes. Due to memory limitations, the 
spatial sizes are chosen to depend on the temperature, the 
largest ones being at the highest temperatures, for which 
Nt is smallest. We find that, in the ghost case, neither 
the propagator nor the dressing function show any visi- 
ble volume or discretization effects within the statistical 
errors. 

The volume dependence of the gluon propagators for 
the two Setups below the thermodynamic transition is 
shown in Fig. 31 As expected, at zero temperature, Dt{p) 
and Dl{p) coincide. In this case no clear volume depen- 
dence is observable, since the three spatial volumes are 
rather similar. At T « 119 MeV, Drip) and Dl{p) are 
already quite different (see discussion below for details). 
On the other hand, these two functions still do not show 
any strong volume dependence. Also note that the con- 
sequences of the violation of rotational symmetry a.re not 
stronger than in the zero-temperature case (23 , [3^ . 



The situation is different when considering the high- 
temperature phase (see Fig.[S]). Indeed, at T « 298 MeV, 
the two propagators show a clear, even though different, 
volume dependence. In particular, the 3d-longitudinal 
gluon propagator is enhanced at small momenta with 
increasing volume, while the 3d-transverse propagator 
becomes weaker in the IR with increasing volume. Fi- 
nally, at T « 597 MeV, the volume dependence of the 
two propagators is similar, i.e. they both become smaller 
(in the IR region) as the volume increases. Moreover, 
the 3d-transverse propagator shows a distinct maximum 
at p « 500 MeV for the two largest volumes, while the 
3d-longitudinal propagator seems to go to a constant at 
small p. These IR behaviors at T w 2Tc are in agree- 
ment with previous studies [3]|. Note that the appar- 
ently stronger finite-size effects seen at T ^ Tc compared 
to the low-temperature cases are likely due to the larger 
spatial volumes (computationally) accessible in the high- 
temperature phase. Also note that a maximum in the 
3d-transverse propagator in the high-temperature phase 
is observed already for a (spatial) lattice volume of 32^. 
At T = 0, in the 4d case, one sees a clear maximum for 
V < 32'* only when simulations are done in the strong- 
coupling regime [13]. Thus, finite-size effects for the 
transverse propagator seem to be smaller in the high- 
temperature phase, in agreement with previous results. 

In Fig. [S] we check for discretization effects by com- 
paring results for three different /3 values for Setups that 
have, within a few percent, the same temperature (i.e. 
T Ri 575 MeV) and the same physical spatial volume 
Vs ~ (3.35 fm)^. The two gluon propagators show a 
different behavior. Indeed, the 3d-transverse propagator 
is only weakly affected by the /3 value at the smallest 
momenta, the effect being of the order of about 20 per- 
cent at most. On the other hand, for the 3d-longitudinal 
propagator, the effects in the IR are more pronounced, 
with variations of up to a factor of almost 2. This in- 
dicates that the scaling limit has not been reached yet 
for these Setups in the IR region. Thus, results for the 
3d-longitudinal gluon propagator at momenta below 1.5 
GeV should be taken with caution. Considering also the 
finite- volume effects discussed above, one probably needs 
to consider larger and finer lattices especially in the case 
of the 3d-longitudinal gluon propagator. 

We can now compare the gluon propagators at differ- 
ent temperatures T. In particular, in Fig. [7|we consider 
Setups with the same physical spatial volume, while in 
Fig. [5] we show data using always the largest spatial vol- 
ume available for a certain temperature. The results are 
similar in the two cases and in the the following we will 
mainly refer to Fig. [S] One clearly sees that at sufficiently 
large momenta the propagators at T > coincide with 
the ones at zero temperature. This is expected, since 
for p ^ T zero-temperature perturbation theory should 
re-emerge. At smaller momenta, the 3d-transverse gluon 
propagator decreases monotonically as the temperature 
increases. Also, there is no sign of sensitivity to the ther- 
modynamic transition. The situation is different for the 
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FIG. 7: The 3d-longitudinal and 3d-transverse gluon propagators and dressing functions as a function of the temperature for 
the same spatial volume 14 of approximately (3.3 fm)"^. We report data for the 3d-longitudinal propagator in the top-left panel 
and for the 3d-longitudinal dressing function in the top-right panel. Results for the 3d-transverse propagator are shown in 
the bottom-left panel and for the 3d-transverse dressing function in the bottom-right panel. Full circles, full triangles, open 
circles, and open triangles correspond to Setups 3 (T = MeV), 6 (T « 119 MeV), 8 {T ^ 298 MeV), and 13 (T « 597 MeV), 
respectively. Momenta respectively below and above 2.4 GeV are measured along the a;-axis and along the spatial diagonal. 
The break at 2.4 GeV is due to violation of rotational invariance. 



3d-longitudinal propagator. Indeed, in the IR region, 
it increases with T when T T^, while it strongly de- 
creases when going from T « to T « 2 T^. From these 
results it is tempting to assume that the 3d-longitudinal 
gluon propagator is largest at the thermodynamic tran- 
sition. However, due to the strong finite-volume and 
discretization effects discussed above, a confirmation of 
this conjecture requires more systematic studies on larger 
and finer lattices. Note that, due to renormalization ef- 
fects, all dressing functions should vanish logarithmically 



for sufficiently large momenta. Clearly this ultraviolet- 
asymptotic regime is not reached yet with our range of 
momenta. 

Finally, from Fig. [5] we see that the Faddeev-Popov- 
ghost propagator does not show any visible temperature 
dependence for all temperatures T considered here, in 
qualitative agreement with previous exploratory studies 
|5l| . We have also studied the spectrum of the Faddeev- 
Popov operator (FPO) as a function of the temper- 
ature T. In particular, we checked if and at what rate 
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FIG. 8: As figure [T] but considering for each temperature T the largest spatial volume available. Full circles, full triangles, 
open circles, and open triangles correspond to Setups 3 (T = MeV), 7 (T ^ 119 MeV), 10 (T ^ 298 MeV) and 15 (T 597 
MeV), respectively. 



the lowest eigenvalue ujs of M vanishes in the infinite- 
volume limit [13]. We recall that, in the 3d case, this 
eigenvalue vanishes faster than the lowest eigenvalue of 
the Laplacian [s^. The behavior of cug as a function of 
the spatial length TV^a = Lg = {VsY^'^ is shown in Fig. 
1101 We fitted our data using a power-law function cx^^ . 
Results of these fits are reported in Table HH We see that 
the exponent h varies little with T and is greater than or 
equal to 2, i.e. the eigenvalue goes to zero at least as fast 
as the lowest eigenvalue of the Laplacian. The magnitude 
c is also essentially unaffected by the temperature. This 
implies that, for the lattice volumes and the lattice spac- 
ings considered here, the eigenvalue ujs (as a function of 
l/Lg) is already in the scaling region. By fitting all data 



we get an exponent h = 2.3(1) (see bottom of Fig. [TO)) . 
to be compared with the exponent of approximately 2.6 
in the 3d case [36j . 

Summarizing our results for the propagators we find 
that 1) the ghost propagator does not have any notice- 
able temperature dependence, 2) the 3d-transverse prop- 
agator depends on the temperature, but it has the same 
(IR-suppressed) behavior below and above the thermody- 
namic transition and 3) the 3d-longitudinal gluon propa- 
gator does show a dependence on the temperature, with 
a different behavior above and below the phase transition 
and an apparent IR divergence at Tc. Moreover, at large 
T, the 3d-transverse propagator is suppressed in the IR 
limit, while the 3d- longitudinal one seems to behave as a 
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FIG. 9: The ghost propagator and dressing function as a func- 
tion of temperature. The top panel shows the ghost propaga- 
tor and the bottom panel the corresponding dressing function. 
Full circles, full triangles, open circles, and open triangles cor- 
respond to Setups 4 (T = MeV), 7 {T ^ 119 MeV), 10 
(T ^ 298 MeV) and 15 (T ^ 597 MeV), respectively. 



massive particle, approaching a plateau for small p. In- 
deed, for all temperatures, volumes and discretizations, 
the 3d-longitudinal giuon dressing function vanishes in 
the IR. This imphes that the 3d-longitudinal propagator 
is always less divergent than that of a massless parti- 
cle. However, in order to check if the behavior is indeed 
that of a massive particle, one needs much larger physi- 
cal volumes. Let us stress that the results obtained here 
at the highest temperature (i.e. T « 2Tc) agree quali- 
tatively with those found in the dimensionally reduced 
theory [sl, [H, [sH] and with previous studies of gluon 
propagators at high temperature [sij . 
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FIG. 10: The dependence of the lowest eigenvalue ujs of the 
FPO on the spatial extension Ls- The lines are fits of the 
type cx~^. In the top figure, we show fits for each temperature 
separately. The fit parameters can be found in Table [III Solid 
lines and circles are at T = MeV, dashed lines and squares 
are at T « 119 MeV, dotted lines and triangles are at T « 298 
MeV while dashed-dotted lines and upside-down triangles are 
at T « 597 MeV. In addition to the Setups listed in Table 
Uwe obtained some of the data using lattice volumes 4x18^ 
and 2 x 22^ at /3 = 2.3. In the bottom figure we consider data 
for all the temperatures. 



A possible explanation of the above results will be dis- 
cussed in Section Ivl 



IV. DYSON-SCHWINGER EQUATIONS 

Due to the strong finite- volume effects observed in the 
numerical data, it is important to complement the pre- 
vious results with ones from the continuum and in the 



13 



TABLE II: Parameters for the fits in Fig. [TO] Fits have been 
done with gnuplot. 



T [MeV] 


c [GeV^] 


h 





0.13(2) 


2.3(2) 


119 


0.119(1) 


2.00(1) 


298 


0.15(1) 


2.31(6) 


597 


0.11(1) 


2.19(5) 



infinite volume, e.g. using DSEs. Let us recall that stud- 
ies at finite T using DSEs have already been presented 
in [H, [13] ■ However, in both cases, in addition to the 
truncations usually employed at zero temperature p^ . 
various approximations have also been performed in the 
solutions of the DSEs. Here, these additional approxi- 
mations will be removed. In this section we discuss the 
(analytic) asymptotic IR solutions of the DSEs at any 



finite temperature, always assuming that external mo- 
menta k obey k <C Aqcd- A (numerical) solution for 
all momenta, using the same truncation scheme as in 
[33I [5^ , is presented in Appendix |D] Let us stress that 
the analysis presented here and in the Appendices IbI - 
ID] applies to any semi-simple gauge group without any 
qualitative change. 

The setup for the system of DSEs at finite T has al- 
ready been extensively discussed elsewhere [H, [H, HI] . 
Let us recall that, in the far infrared, one can use the so- 
called ghost-loop-only truncation [53], which keeps only 
the leading infrared term arising in the derivation of the 
DSEs from the Faddeev-Popov-ghost term in the gauge- 
fixed Lagrangian. In the Gribov-Zwanziger scenario, this 
term dominates the action at low momenta. 

Then, the equations for the zero (soft, fco ~ 0) modes 
are given by [5 31 ] 



G(0,/c) 
1 

Z(0,fc) 
1 

7J(0, k) 



(27r)3 



^ / dqde \AT{0,qo,k,q)G{qo,q)Z{qo,q-k) + Al{0, qo,k, q)G{qo, q)H{qo, q- k)] (10) 



^3T 



9'TCa 



3L 



/ dqd0RiO,qo,k,q)G{qo,q)G{qa,q + k) 
^' qo •' 

^^T.J dqd0P{O,qo,k,q)G{qo,q)G{qo,q + k) . 



i2n) 



(11) 
(12) 



Here, we use the definitions G(fc^) = k'^Daik'^), Z{k'^) ^ 
k^Drik^), H{k^) = k'^DLik^), where Daik^), Drik^) 
and DLik"^) are, respectively, the ghost propagator and 
the transverse and longitudinal gluon propagators. The 
integral kernels At, Al, R, P and the angle 6 are defined 
in Appendix iBl Also, T is the temperature, g is the bare 
coupling constant and Ca indicates the adjoint Casimir 
of the gauge group [for SU{2), Ca — 2]. The graphical 
representation of this system of equations is given in Fig. 
[TT] Note that Eqs. (jlllfT^ have been obtained from the 
tensor equation for the gluon propagator by contraction 
with generalizations of the projectors defined in Eqs. ^ 
and ([3]). These tensors are parameterized by the variables 
C and ^ |33j , which appear explicitly only in the integra- 
tion kernels R and P (see Appendix |B| . Variations of 
these parameters can be used to investigate truncation 
artifacts and, in particular, the appearance of spurious 
divergences.^ Let us also note that the equations above 
reduce to the corresponding equations of a 3-dimensional 
Yang-Mills-Higgs system when the hard modes {qo ^ 0) 
inside loops are neglected [l^. Then, due to the ghost- 



loop-only approximation, the Higgs field H{0, k) decou- 
ples'' and G{0, k) and Z{{), k) become the dressing func- 
tion of a pure 3d Yang-Mills theory [l^]. In Appendix 
[Cl we show that the solutions of Eqs. pUlfT^ at zero 
temperature coincide with the well-known solutions pre- 
sented in Refs. [l^. Finally, we remark that one can find 
different wave-function renormalization constants {Z^t 
and ^3^) in the gluon equations due to possible finite 
contributions arising at non-zero temperature. 
Equations pTUFT^ can also be written as 



— i-^ = ^3 + HG(0,fc) + V HG(go,fc) (13) 
G(0, k) 



90 #0 



— i-^ = Z^T + Hz(0,^) + 112(90, fc) (14) 
Z(0, k) 



go#0 



(15) 



where we indicate with H^ the various self-energies. Note 
that, in the 3d- longitudinal equation the zero-component 



See e.g. [Tsl. ISSl lssIl for a detailed discussion of this topic at finite 
temperature. 



^ Indeed, it is easy to check that for fco = go = one finds Ai^ = 
P = 0. 
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FIG. 11: The DSEs in the ghost-loop-only approximation at 
finite temperature. Dotted lines represent ghosts, curly lines 
are used for the 3d-transverse gluons and dashed lines repre- 
sent 3d-longitudinal gluons. Lines with a dot indicate a full 
propagator instead of a bare one. Vertices with a small dot 
are bare ones, while those with a large, empty dot represent 
full vertices. However, in our calculations these vertices are 
also taken to be bare. 



more, the ghost-gluon vertex is taken equal to the bare 
one. This approximation is supported by various calcu- 
lations at zero [H, [H, [l^, [s^l and at infinite temperature 
(25l[36j. Of course, when considering the far IR limit at 
non-zero temperature one needs to ensure the condition 
<C T for the external momenta. In this limit one ob- 
tains [3^ that the hard-mode dressing functions Z{ko, k), 
H(ko,k), G(ko,k) become constant in the infrared, i.e. 
they have the behavior of the dressing function of mas- 
sive particles. Thus, for fc ^ 0, the dressing functions 
of the hard modes for the 3d-transverse gluon, the 3d- 
longitudinal gluon and the ghost, are respectively given 
by ^2(^0), Ah{ko) and Ag{ko), which depend only on the 
Matsubara mode feg. In order to find a solution to the 
system of equations (ITU] [T2), we need these IR constants 
to be bounded. 



Finally, for the zero-modes (fco 
power-law ansatze 



0) we consider the 



qo — explicitly vanishes, since P{0, 0,k,q) — 0. This is 
a direct consequence of the tensor-structure of the ghost- 
gluon vertex. Thus, the behavior of the longitudinal 
mode H{0, k) is entirely determined by the hard modes 
and depends only implicitly on the soft ones. 

The treatment of the full system, given in Appendix 
IdI justifies a-posteriori this ghost-loop-only truncation 
for momenta k in the far IR limit. The same trunca- 
tion is also sensible at zero temperature [13] • Further- 



G(0,fc) = Bgk"^" 
Z{0,k) ^ B^k^^ 
H{0,k) - Bhk^' 



(16) 
(17) 
(18) 



which depends on the IR exponents k, t and I and on the 
constant coefficients Bg, Bz and B^. Then, for all values 

of K and t, the self-energies 110(0, fc) and 11^(0, fc) take 
exactly the same form as in three dimensions [ISj], i.e. 



nG(0,fc) = -BgBzg^TCA 



T(2 + 2k) 



K(3 + 4(-2 + K)K)r (2k-|-|)' 



,,t + K- 



n (n A r2 2-4(i+«) (-2 + C + 2>,{C - l))r(2 + 2k) sec(2^Ac) sin(^A.)^ 

Tlz(0,k) = B g TCa . . , 3. V 

K^(l -I- k)1 (2k 4- I) 



(19) 
(20) 



where y = k^ = k^. Clearly, if these contributions are not exactly canceled (or dominated) by the remaining Matsubara 
sums, they will give the leading IR behavior of G(0, k) and of Z{0, k). In this case one finds a 3d-type behavior for all 
non-zero temperatures T. To show that this is indeed the case we observe that — by setting the hard-mode dressing 
functions equal to the constants Ag{ko), ^^(fco) and Ah{ko) — the self-energies in the ghost equation (|13p can be 
rewritten as 



IlG{qo^O,k) = g^TCA 



Aziqa) , Ahiqo) 



12n 



96 



1 

M 



Az{qa) , Ah{qa)\ k 



SOtt 



1920 



l9o| 



+ k^TT{qo,k) 



(21) 



Here, the first term is logarithmically divergent and it 
can be absorbed in the renormalization constant. At the 
same time, the second term is sub-dominant when com- 
pared to the 3d-term, for values of k and t permitted by 
integral convergence [isj. and is finite after summation 
over (Jo- Finally, the function n{k) vanishes identically as 



k ^ 0. Thus, the leading part of the ghost equation is the 
same as in the 3d-case, i.e. it is given by nG(0, k), pro- 
vided that the quantities Az{qo) and Ah{qo) do not rise 
strongly as qo goes to infinity. Actually, due to asymp- 
totic freedom, these quantities vanish logarithmically in 
the ultraviolet limit. 
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Finally, we should check that a 3d-type behavior is ob- 
tained for all non-zero temperatures T also when consid- 
ering the gluon equations pT|l and (|12p . To this end, we 
first discard spurious divergences in these two equations. 
These are quadratic divergences, vanishing for ^ = 3 in 
the 3d-transverse equation ()11|) and similarly in the lon- 



gitudinal equation ((T2l) . One can show that they are ar- 
tifacts of the truncation scheme considered [l^, [H, [11] . 
After subtraction of these divergences, one can write the 
contributions of the hard-modes in the self-energies [see 
Eqs. (HlHlSl)] as 



T^z{qa,k) = , n^((7o,fc) ^ ^ . (22) 

In both cases the fmiction tTj vanishes as /c ^ 0, separately for each Matsubara term. Thus, the subtracted part 
cannot contribute to the IR behavior in the 3d-transverse equation. This is the case* also after summing over qq . On 
the other hand, the IR contributions from the un-subtracted self-energies cannot in general be neglected, as done when 
considering only the contributions (|^. This requires, of course, a regularization [1^ of the spurious divergences.^ To 
achieve this, one can replace the approximately constant dressing functions of the hard modes by Ag{qo){q'^ + qo)~^j 
which are clearly suppressed when r > 0. This is the prescription commonly used for regularizing the divergences at 
zero temperature with r = —k pTllTsl [l6l[T9j . Then, the integrals can be performed and in the limit fc ^ one finds 

nf(.o,A:^0) = -^|^|^|-^(2.T)--^^,(n)>|-- (23) 

^ ' n#0 



where n = qo/2'KT. These sums diverge for r < 1/2, 
due to the term (?q~'*^. For r = 1/2, this exponent be- 
comes equal to —1, i.e. the sums are logarithmically di- 
vergent. Finally, for r > 1/2, the sums are finite and 
they can (in principle) be resummed analytically. For 
example, for Ag{n) — (27rT)^'^ one finds that the sum 
m n IS equal to 2(27rT)''^Ci? (4r - 1), where Cr is the 
Riemann-C function. One should also note that the terms 
in (P5HM)) behave like l/k'^, i.e. like a mass term in the 
IR limit. In the 3d-transverse case — but not in the 
3d-longitudinal case — this term may not be renor- 
malized, as this is not allowed by gauge-invariance. If 
one does not allow either a divergence stronger than log- 
arithmic (even though this cannot be excluded by the 
perturbative renormalization theorems) or different ex- 
ponents T to the 3d-transverse and to the 3d-longitudinal 
case, then the only possibility is to have r > 1/2. For 
all such values of r, the contribution (|23p is sub-leading 
in the 3d-transverse equation (HH). At the same time, 
the screening mass in the 3d-longitudinal equation would 
then be solely due to the regularized contribution ((24)) . 
This result is a consequence of the truncation scheme. 
In fact, due to the decoupling of the soft modes in the 
3d- longitudinal equation [see Eq. (fT5|) ]. the result is dom- 

* Note that this also applies to the remaining gluon loops, which 
are not treated explicitly here. 

® Note that, in the 3d-transverse case, divergences appear for each 
hard mode, while in the 3d-longitudinal case only the sum over 
the hard modes is affected by spurious divergences. 



inated by the hard modes, which are very sensitive to 
truncation artifacts, since they live on a scale that is ef- 
fectively mid-momentum. Thus, the present truncation 
scheme is not able to yield a consistent description of 
the electric screening mass and a determination of its 
value is not possible. Nonetheless, the 3d-longitudinal 
gluon propagator seems to exhibit a screening mass, and 
thus gives likely a qualitatively correct description of the 
physics involved. This screening mass would not allow 
the 3d-longitudinal gluon propagator to modify the IR 
behavior of either the 3d-transverse gluon propagator 
or of the ghost propagator. Therefore, as shown above, 
these two propagators should behave at all non-zero tem- 
peratures exactly as in the three-dimensional case when 
one considers momenta much smaller than the temper- 
ature T and than Aqcd- It is then plausible that for 
momenta in the range T ^ p <C Aqcd the gluon and the 
ghost propagator would still exhibit a behavior similar 
to that found in the four-dimensional case. This is also 
found in DSE studies in a finite volume [isj . 

The behaviors described in this section are in agree- 
ment with calculations using renormalization-group tech- 
niques in a background- field gauge [s^ . Indeed, when 
going from T = to T > 0, also in that case one 
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finds a discontinuous change in the behavior of the run- 
ning coupUng, i.e. in the IR region the running cou- 
phng switches from a four-dimensional to an effectively 
three-dimensional behavior. Let us also note that previ- 
ous studies using DSEs [52] had assumed that the zero- 
temperature behavior persists up to the thermodynamic 
transition. Here we have relaxed this hypothesis and 
eliminated some of the approximations of the numer- 
ical method employed in }52l |. showing that the high- 
temperature results apply also to non-zero temperatures 
below the critical temperature T^, a phenomenon previ- 
ously interpreted as super-cooling [33| . 

V. SUMMARY 

Let us summarize here the results obtained above for 
the 4d finite-temperature case. From our lattice calcula- 
tions we have that: 

• The 3d-transverse gluon propagator decreases as 
the temperature increases. It also seems to have 
smaller finite-volume effects (in the IR limit) and 
stronger IR suppression at high temperature than 
at zero temperature. 

• The 3d-longitudinal gluon propagator at small mo- 
menta increases as the temperature goes from to 
Tc, seems to diverge at Tc and then drops as T 
becomes larger, apparently reaching a plateau for 
small p. Also, in the IR limit, it is less divergent 
than the propagator of a massless particle, with a 
possible exception near the phase transition. 

• The ghost propagator is nearly temperature-inde- 
pendent. 

• The smallest eigenvalue of the Faddeev-Popov op- 
erator, considered as a function of the spatial lattice 
side Ls, goes to zero faster than 1/L^ as Lg goes to 
infinity, for all temperatures. 

From the DSEs we see that: 

• The IR behavior of the gluon and ghost propaga- 
tors changes abruptly from zero to any non-zero 
temperature. 

• For momenta p <C (T, Aqcd), the ghost and the 
3d-transverse gluon propagators show the same be- 
havior obtained in the dimensionally reduced the- 
ory [ll,[33| and the IR behavior in the spatial sector 
is in accordance with the prediction of the Gribov- 
Zwanziger scenario [HI, [13] . 

• The 3d-longitudinal gluon is likely to have a dy- 
namical mass at any non-zero temperature. 

• These results can be connected continuously to per- 
turbation theory at large p (see Appendix |D]) and 
to the solution at zero temperature (see Appendix 

El). 



As explained in Section IIVI above, the main assump- 
tion for these DSE results is that the electric screening 
mass is non-zero at all temperatures. Indeed, in the case 
of a null screening mass, one can find different solutions 
for the DSEs. In particular, the ghost propagator may or 
may not be IR divergent, while the 3d-transverse gluon 
propagator would still be IR suppressed [53*]. Note that 
the analytic asymptotic results from the DSE calcula- 
tions are valid for momenta p <ti min(T, Aqcd)- As al- 
ready stressed in Section IIII A| this range of momenta 
is not accessible with the lattice volumes used here. At 
intermediate momenta, DSE results agree qualitatively 
with the lattice results for all temperatures considered. 

We can try to interpret these results by considering 
the Gribov-Zwanziger scenario applied to the 4d T = 
case and to the dimensionally reduced theory. As ex- 
plained in Section IIII A[ this scenario is supported by 
lattice [p, [H, Hi, [si, [H, HI] and by DSE calculations 
[Til [l3l.Tl9|. Then, at T > 0, due to the compactifica- 
tion of the time direction, the configuration space near 
the Gribov horizon should be essentially reduced to that 
of a three-(space-tinie)-dimensional system and, in the 
infinite- volume limit, the configurations should lie on the 
Gribov horizon [13] for all temperatures and the Gribov- 
Zwanziger confinement scenario would apply to the so- 
called magnetic sector. ^'^ As a consequence, one can re- 
move the so-called IR problem [i, [3, [H, [11] since an in- 
frared suppressed 3d-transverse gluon propagator cancels 
the pcrturbative infrared divergences. 

Based on the above considerations, our results may be 
organized into the following scenario: 

• At non-zero temperature, the time-like momenta 
are no longer continuous, but discrete. Due to 
this gap, the properties of configurations near 
the Gribov-horizon are changed. This modifies 
(probably only quantitatively) the spectrum of the 
Faddeev-Popov operator. Thus, the enhancement 
for small eigenvalues observed at T = [s^ should 
still be present and the IR enhancement of the 
ghost propagator might be reduced but not elim- 
inated. This is verified by our lattice data. 

• The consequences for the 3d-transverse gluon and 
the 3d-longitudinal gluon are different, due to the 
vector character of the ghost-gluon interaction. In 
particular, in the 3d-trans verse case, the coupling 
is not modified in the DSEs and this propagator 
is still IR suppressed. Moreover, this suppression 
is actually stronger due to the structure of the 3d- 



This is the case not only in the Landau gauge but also in Coulomb 
gauge l57l Issll . Note that this is in line with the fact that, in 
the spatial subspace, the interactions are qualitatively unaltered 
when the temperature is switched on and all purely spatial ob- 
servables, such as the spatial string tension 8], are not modified 
by the transition. 
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transverse gluon-ghost interaction. This is also sug- 
gested from our lattice data. 

• On the other hand, the coupling of the ghost to 
the 3d- longitudinal gluon becomes gapped. Thus, 
the solution of the DSEs is consistent with these 
gluons acquiring a real (dynamical) mass, being no 
longer IR suppressed. This 3d-longitudinal screen- 
ing mass is generated solely by the hard modes, 
leading likely to a mass of the order of the tem- 
perature, which is the characteristic scale of the 
hard- modes. In addition, this is observed in lattice 
calculations at very high temperatures [Slj . At the 
same time, this screening mass should, similarly to 
the hard modes, be sensitive to the phase transi- 
tion. This is partially seen from our lattice data, 
since the longitudinal propagator is sensitive to the 
phase transition and seems to reach a plateau at 
small p for T > Tc. 

• Finally, it should be noted that the presence of an 
electric screening mass does not imply that the 3d- 
longitudinal gluon is deconfined or an observable 
particle. Indeed, a massive particle does not nec- 
essarily have a positive-definite spectral function 
[13L [33, 153] ■ Moreover, the Ao-field is (at the per- 
turbative level) member of a BRST-quartet at any 
finite temperature [60| . These considerations would 
imply that no gluon is part of the physical spectrum 
also at any T > 0." 

A number of predictions follows from the scenario pre- 
sented above, which could in principle be further tested 
using lattice calculations at sufficiently large volumes and 
with sufficiently fine discretization. The main predictions 
are: 

• For all temperatures, the 3d-transverse gluon prop- 
agator is IR suppressed, while the ghost propagator 
is IR enhanced. 

• When p < (TjAqcd) the IR behavior of the 3d- 
transverse gluon propagator and of the ghost prop- 
agator is quantitatively different from the behavior 
observed for T <C p ^ Aqcd, the latter one being 
similar to the behaviors observed at zero tempera- 
ture with momenta p <ti Aqcd ■ 

• The 3d-longitudinal gluon propagator is null at zero 
momentum only at T = 0, while it goes to a fi- 
nite (non-zero) value a,t p — for all temperatures 
T > 0. This propagator is sensitive to the phase 
transition. 



Of course, it would be interesting to compare our nu- 
merical results to a similar study for the SU{3) group, 
for which a different kind of thermodynamic transition 
is expected Note that first investigations, with and 
without dynamical quarks, have been performed in (63i |. 

VI. CONCLUSIONS 

We have presented an analysis of gluon and ghost 
propagators at finite temperature, using lattice gauge 
theory and DSEs. Results using these two approaches 
seem to be (at least at the qualitative level) consistent 
with each other. In particular, when the temperature is 
turned on, one finds different effects for 1) the tempo- 
ral sector, which includes the soft 3d-longitudinal gluon, 
2) the spatial sector, containing the soft ghost and the 
soft 3d-transverse gluon and 3) the hard sector with the 
(essentially perturbative) hard modes. In the spatial 
sub-sector, the IR behavior is similar to that of a 3- 
dimensional theory and the correlation functions agree 
with those expected from the Gribov-Zwanziger confine- 
ment scenario. Also, their dependence on the temper- 
ature is smooth. On the contrary, the temporal sec- 
tor is sensitive to the thermodynamic transition, being 
screened and decoupling at the scale of the temperature. 
Nevertheless, our results are consistent with the possibil- 
ity that all gluons are confined at all temperatures. Note 
that this does not alter the bulk thermodynamic prop- 
erties, which still have a Stefan-Boltzmann-like behavior 
[H, m, [13, ■ These results can be understood by con- 
sidering the Gribov-Zwanziger confinement scenario at 
T 7^ 0. As a consequence, gluon confinement is not af- 
fected by the thermodynamic transition. This would ex- 
plain the presence of non-perturbative effects in the high- 
temperature phase and the confining properties of the 
dimensionally-reduced theory in the infinite-temperature 
limit. 

One should, however, recall that these results are af- 
fected by various (technical) limitations. Indeed, much 
larger and finer lattices and more sophisticated DSE 
schemes are needed in order to obtain a complete under- 
standing of the effect of the temperature in the gluonic 
sector of Yang-Mills theory. 

Finally, let us note that the introduction of dynamical 
quarks should not modify the Gribov-Zwanziger scenario 
and, as a consequence, the scenario described above. At 
the same time, the chiral transition should also not be 
affected by the IR arguments considered here, as the rel- 
evant dynamics occurs at the scale of Aqcd (see e.g. [T6|). 
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APPENDIX A: ASYMMETRIC LATTICES 

When using asymmetric lattices at T = (e.g. if 
the time extension of the lattice is larger than the spa- 
tial ones) it is convenient to consider the two scalar 
gluon propagators given by D±{p) = Dii{p) / {d — 1) and 
D\\{p) = Dooip), where in the first case we sum over the 
spatial indices i only. In Fig. [12] we report the results ob- 
tained for these two different tensor components of the 
gluon propagator when considering an asymmetric lat- 
tice of size 24^ x 128 (at /3 = 4.2). Let us note that 
the total volume in this case is roughly the same as for 
a 40'^ lattice. In this symmetric case, one clearly sees a 
maximum in the propagator for a momentum p « 400 



MeV (see Section UlL^ and the top plot in Fig. [J). In 
the asymmetric case, on the contrary, it is evident that 
the two tensor components show a very different behav- 
ior and that the behavior of the orthogonal propagator 
D±{p) depends on the type of momentum considered, i.e. 
along the short axis of the lattice (ps) or along the long 
axis of the lattice {pi). Also, there is no clear evidence of 
a maximum, even when considering D±{pi). Moreover, 
as observed already in [4l|, the IR suppression obtained 
using an asymmetric lattice N^~^Nt with Nt ^ Ng is 
usually smaller than the suppression found on a symmet- 
ric lattice with volume Nf. This effect clearly reduces the 
advantages of using asymmetric lattices. 

When considering asymmetric lattices, one might also 
use asymmetric couplings in order to have the same phys- 
ical extent in all directions. This would probably reduce 
the systematic effects observed in [4l| , while keeping the 
numerical costs of the simulations small. 



APPENDIX B: INTEGRAL KERNELS 



Here we define the integral kernels obtained with the ghost-loop-only truncation [see Eqs. p0|) ~ p2)) ]. We use the 
following notation 



y = P = kl + P; x^q^ = ql + q^\ k=\k\; q ^ \q\ 
u = X + y + 2fco(7o — '2k ■ q; v = {k — q)^ 
z = X + y + 2koqo + 2k ■ q; w = {k + cf)^ 
9 = arccos ( k ■ q/kq ) . 



Then, the kernels in the ghost self-energy pI7|) are given by 



koq^ + k^qo - (ko + qo)k ■ q 



sm ( 



xyvv? 



At the same time, in the 3d-transverse gluon equation (jlip . the kernel is equal to 



(Bl) 



AT{ko,qo,k,q) = ^^—i , (B2) 

xyuv 

^ 1 2 



AL{ko,qo,k,q) = ^ . (B3) 



g3 [-g + /c(C- l)cos6l-f gCcos6l2] sin6' 

R[ko,qo,k,q) = ^ , (B4) 

Zxyz 



while in the 3d-longitudinal equation (|12p it is 

P{ko,qo,k,q) = -L-[q\qo{-2kliko + qO)iC - 1) + k^[ko - ^{ko ^ qo)]} 

-kqko{ko{C - I) + 2qa{2^ - I)} cos9 + klq^^cose^)]sin0 . (B5) 

Note that the 3d-transverse equation (fTT|) has been contracted with a tensor parameterized by a variable ^ [1^ . This 
variable appears only inside the integration kernel (|B4[) . The case C = 1 corresponds to the projector given in Eq. 
0. The same has also been done for the 3d- longitudinal equation (fT^ . with the variable Again, it only appears 
inside the integration kernel (jB5p . For the zero mode (/cq = 0), ^ factorizes in all terms and can thus be eliminated 
from the equation. Hence the equation depends only implicitly on ^, due to the hard modes [33| . 
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FIG. 12: The orthogonal D±(p) and the parallel -Dy (p) components of the gluon propagator evaluated using the (three- 
dimensional) strongly-asymmetric lattice 24^ x 128 at P — 4.2. Open symbols correspond to momenta (ps) along one of the 
short sides of the lattice, while full symbols correspond to momenta (pi) along the long side of the lattice. Here we considered 
420 configurations and the same settings used for Setup 15. Note that, due to the Landau gauge condition, one has {pi, 0) = 
for all momenta pi . 



APPENDIX C: VACUUM SOLUTION 



As said in Section lWl above. a necessary requirement is that Eqs. (|10H12p reproduce the well-known vacuum behavior 
[ipl at T = 0. To show this, it is useful to consider linear combinations of Eqs. pip and (fT^ . In particular, we consider 
a combination with weights given by 2/3 and 1/3. We also introduce a second combination with weights equal to 1 
and (C — 3)/2. At T = the Matsubara sums become integrals. After performing these integrals in the equations 
obtained by considering both linear combinations, we find 



_ 2 1 _ 2 o2 4-^-^"^^(l + ^)(3 + 2k){-6 + C(3 + 4ac)) csc(27r^) 

3^ 3iJ s 3r(i-K)2r(K + 



Uh^ lC-31 „ pik) 



with y = . The function p{k) vanishes for fc — s- 0. 
Hence, for k > 0.5, l{k) is sub-leading compared to /(fc). 
This yields 



Z(fc) 
Hik) 



2(C~4) 



3(C-3)/(fc)-2Z(fc) 

4-C 
3/(fc)-2?(fc) ■ 



(C3) 
(C4) 



Without considering the trivial case C = 4 and neglect- 
ing the sub-leading quantity l{k), these two functions co- 
incide only in the case C = 1- For other values of ( 
the dressing functions Z{k) and H{k) will have the same 
leading IR behavior, but they will differ by a constant fac- 
tor. This is simply an artifact of the non-0(4)-invariant 
projection considered. Thus, for the remaining discussion 
we can set C equal to 1. 



By using the above results and the IR ansatz (fTB]) , the 
ghost equation gives the relation 



37rr(-2K)r(2 + n) csc{ttk) 
k2(2 - 3k + K2)r(_K)3r(2 + 2k)- 



(C5) 

Here and henceforth the sub-leading terms proportional 
to p{k) have been dropped. Clearly, this relation can be 
satisfied for all values of y (in the infrared limit) only if 
the four-dimensional relation 119, 1201 



2k + < = 



(C6) 



is satisfied. By dividing equations (jCip and (jC5|) . insert- 
ing the power-law ansatze P^HTSl) and setting Z = H , 
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FIG. 13: The truncated Dyson- Schwinger equations at finite temperature. The wiggly lines represent transverse gluons, the 
dashed lines longitudinal gluons and the dotted lines denote ghosts. Lines with a full dot represent self-consistent propagators 
while small dots indicate bare vertices. The open-circled vertices are full and must be constructed in a given truncation scheme. 
A bare ghost-ghion vertex and modified bare gluon vertices have been used here. 



one finds a consistency condition for k, i.e. 

(2-k)(k- 1) 



1 = 



12(3 + 4k(2 + k)) 



(C7) 



This equation is identical to the T = equation and is 
solved by K « -0.595353 [ill,!!!. 



APPENDIX D: THE FULL SYSTEM 

In this section we consider the system of DSEs trun- 
cated at the one- loop level (see Fig. [13]). This truncation 
has already been discussed in detail in Refs. [H, for 
the finite (i.e. nonzero) temperature case. 



1. Perturbation theory and truncation artifacts 

Many direct experimental verifications of QCD are per- 
formed in a regime where perturbation theory is appli- 
cable. Thus, any non-perturbative treatment must make 
contact with it. In the case of DSEs (at the presented 
truncation level) this implies recovering the resummed 
leading order perturbation theory [l^,[la,[i^. Of course, 
at finite temperature T, the vacuum perturbation theory 
is expected to be recovered only at momenta large com- 
pared to T. 

In order to show that this is the case, the first step is to 
identify and remove the spurious quadratic divergences 
[Tsl [To . [Tgj . Since the summation over the Matsubara fre- 
quencies extends to infinity, spurious divergences might 
depend on the order of integration and summation. Fur- 
thermore, it is possible that a contribution is finite with 
respect to integration but not to summation and vice- 
versa. Of course, only the region corresponding to large 
n is relevant for the analysis of spurious divergences. In 



this limit, the difference between two consecutive Mat- 
subara frequencies becomes negligible and the summa- 
tion becomes equivalent to an integration. One can then 
replace | (fl by g sin ij and qo by q cos 7] (with q^ — q^^ (f') 
in the integral kernels K{q) appearing in DSEs. Since 
the angular integration over 77 yields no divergences, it 
is sufficient to replace all integral kernels K{q) by the 
prescription 



K{q) ^ K{q) - 1 hm q^K{q) 



(Dl) 



The resulting integrals have no spurious quadratic diver- 
gences, but only the usual logarithmic one. One could ab- 
sorb the subtracted part in gauge non-invariant counter- 
terms, as they originate purely from violation of gauge 
invariance. On the other hand, as discussed in Section llVl 
above, this affects the IR behavior of the solutions. To 
surpass this problem constructively, the electric screen- 
ing mass is here explicitly renormalized to a fixed value 
nir = g^T, since this is allowed by gauge invariance [ssj . 
Note that this prescription is sub-leading in the ultra- 
violet and does not affect perturbation theory. At the 
same time, the corresponding terms in the 3d-transverse 
equation can be dropped, as they are sub-leading both in 
the IR and in the ultraviolet limit. The only exception 
is the soft-mode contribution, which is treated as in the 
infinite-temperature limit [l3|. As for the genuine tad- 
poles, they could in principle give a finite contribution. 
However, when removing the quadratic divergences ac- 
cording to (|Dip . they vanish identically. Thus, also at 
finite temperature, the tadpoles contribute only a pure 
divergence, which can be absorbed in gauge non-invariant 
counter-terms. Also note that this implies that all spuri- 
ous divergences are not affected by a finite temperature.^^ 



This could be the case, for an inconsistent truncation scheme. 
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This is to some extent not surprising, as the usual diver- 
gences in quantum field theory may not be affected by a 
finite temperature (60| . 

As a second step, we have now to fix the renormal- 
ization prescription for the (physical) logarithmic diver- 
gences. In order to recover the 3d-theory in the infinite- 
temperature limit, two aspects have to be considered in 
the renormalization prescription ^33j . Firstly, one needs a 
finite 3d-coupling. To this end, one requires the quantity 
g^T to be fixed to Aqcd for large temperatures [33l.[53l|.^^ 
This implicitly defines a temperature-dependent renor- 
malization scheme and automatically fixes the renormal- 
ization scale /z as a function of T by the running of the 
coupHng constant g^/{^iT) = as- At the same time, a 
smooth zero-temperature limit is obtained using the pre- 
scription 



M = Aqcd { exp 



/3o 



1 -I- exp 



^oAqcd . 



(D2) 

Here /3o is the first coefficient of the /3-function. Then, 
according to the perturbatively resummed 1-loop result, 
the largest value attained by is 1, i.e. = An. 

Secondly, the wave function must be renormalized such 
that the dressing functions become unity at infinite mo- 
mentum for T — > GO. This is guaranteed by selecting the 
subtraction point 



= lOU 



QCD 



27rT 



(D3) 



Note that the term W^Aqqd has been added to renor- 
malize perturbatively at zero temperature. Clearly, the 
choice of temperature-dependent values for s and jj, en- 
tails that the renormalization constants will also contain 
finite temperature-dependent contributions. 

The regularization and renormalization are performed 
as in Refs. [l^jllll by adding explicit counter-terms. This 
makes multiplicative renormalizability manifest. The 
renormalization conditions are 

G{s) = Z{s) = H{s) = 1 . (D4) 

'^'^ The precise value of A^qcd is irrelevant for the present purpose 
but, as for the zero temperature case [Tsll . it could be obtained 
from comparison with the running coupling at sufficiently large 



They implement the requirement 

G'^{s)Z{s) = G^{s)H{s) = 1 



(D5) 



from Slavnov- Taylor identities (STI) at zero temperature 
im . Note that this prescription does not allow any free- 
dom in the values of the IR pre-factors in the ansatze 
(ll6H18p . This renormalization scheme is explicitly inde- 
pendent of any Matsubara frequencies, and thus it suf- 
fices to determine the wave function renormalization con- 
stants in the soft equations. On the other hand, this 
renormalization scheme only works if performed at mo- 
menta sufficiently far in the perturbative regime, since g 
has only been included to the leading perturbative order. 

In order to recover perturbation theory, the ultra- 
violet asymptotic dressing functions must be given by 
the renormalization-group-improved perturbative propa- 
gators t66|] 



G{P) 


= Gis) 


wlog 




+ 1 


H{p) 


= Zis) 


UJ log 




+ 1 



4(7+1) 
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(D6) 

(D7) 

(D8) 
(D9) 



This also implies that the longitudinal and transverse 
gluon propagators must coincide at sufficiently large mo- 
menta p. At the same time, any kind of gauge-invariance 
violation must vanish, i.e. the propagators should become 
independent of the variables C and ^. 

It thus remains to test whether the truncation scheme 
permits this. At T = it has been found that only a 
dressed three-gluon vertex yields the correct perturba- 
tive solution [19|. Thus, we consider the Bose-symmetric 
ansatz 

(perturbative) momenta. 



r(p,g,%,p - A{q,p,k)T{p,q,k)fJp{p,q,k) (DIO) 
= aG(^{q^+p^ + k^)^ \(^{q^+p^ + e)^ \{p,q,k)fX{p,q,k) . (Dll) 

Here, T and T^") are respectively the full and the tree-level vertex and the constants a, ac and az will be chosen 
such that the propagators (|D6HD7[) are a solution of the system. Of course the above ansatz is valid only at high 
momenta. The low- momentum case, for which Z ^ H , will be discussed below. 



As usual, due to the STI and the results obtained in the non- 
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Finally, it is sufRcient to inspect the limit p ^ po in a, given equation, which thus reduces to the corresponding 
soft equation. As only the large-momentum behavior is relevant, the integrals can be treated with a lower cutoff p. 
Furthermore, the denominators can then be expanded in inverse powers of the integration momentum^^ up to 0{q^^). 
Approximating further internal dressing functions by G{p + q) fa G{q) etc. p^ . a simple set of equations is obtained 

1 ,B.) 



G(p2) 647r2 " y 

- ^3T + / dy— / dy (D13) 



Z(p2) 3847r2 " y SSAtt^ Jp2 " y 

C 9^Ca^ [""^ G{yf 25g2c^^ A{y,y,p^)Z{yf 

- C^3L - -TT^ / dy— ——^ / dy . (D14) 



(The integral kernels and other details of these equations 
can be found in [ssj].) In these equations, the asymptotic 
form (jD7|) has already been used to replace H hy Z and 



y = q^ = q^ + (f' ■ Note that the third equation is the 
original 3d-longitudinal equation and the ^-dependence 
can obviously be removed. The wave function renormal- 
ization constants are chosen to cancel exactly the upper 
boundary of the integrals, thus rendering the equations 
finite and independent of the regulator SI. 

One can then check that the ghost equation (|D12[) is 
solved by (|D6HD7p . This yields the relation (jD9p between 
7 and 5 and the value of to [see Eq. (|D8|) ] . On the other 
hand, the value of 5 is not fixed. The remaining task is 
the choice for A{y,y,p^). Using the ansatz (jDll[) in the 
3d-transverse equation (|D13[) . one finds for the parame- 
ters az and a the conditions 



az 



-2~6S + aGS 

1 + 25 
-18-41(5 + 3(5C 
5(47 + 30 



(D15) 
(D16) 



At the same time, the equation becomes independent of 
aa- In the 3d-longitudinal equation (|D14[) one finds again 
for az the solution (jPlSp . On the other hand, a is deter- 
mined to be given by 



-9-195 
256 



(D17) 



Incidentally, there is only one value of 5 that at the 
same time removes the Odependence in the 3d-transverse 
equation and yields the same a value in both the 3d- 
transverse equation and in the 3d-longitudinal equation. 
This value is (5 = —9/44, which coincides with the re- 
sult from the renormalization-group-improved perturba- 
tion theory (|D9[) . For this value, one finds a = 1 in both 
equations above. Thus, the condition of gauge invariance 
uniquely requires for this truncation scheme the correct 

perturbative regime [261 . the ghost-gluon vertex is kept bare. 
Originally, at T = 0, the angular integrals were solved exactly 
instead of considering this expansion 11911 . However, due to the 
non-Euclidean invariant projection, this is not as simple any- 



value of 5. This is quite convenient, although likely acci- 
dental. 

Let us note that it is always possible to rewrite a = Zib 
with an arbitrary value for Zi, thus realizing the STI 
Zi = Z3/Z3 deliberately by fixing b appropriately. 

Finally, as said above, we can consider the low- 
momentum form of the three-gluon vertex. To this 
end, we replace G°^Z°^ in Eq. ljDll[) with the quantity 
QnaGT/3+maGL/3znaz/3H^naH/3^ jjcrc, u and m are the 

number of 3d-transverse and 3d-longitudinal legs. The 
exponents az and an are chosen to satisfy (|D15p for 
og = O'er and ac — aQi^, respectively. A convenient 
choice for the remaining parameters uqt and acL is 



acT 



2t{l + 35) 
t5 + K + 25k, 
21(1 + i5) 



15 



25k 



(D18) 
(D19) 



where k, t, and / are the IR exponents of the dressing 
functions PWTS)) . This yields an IR (positive) constant 
three-gluon vertex. Although this is not in agreement 
with recent results on the IR behavior of this vertex [2^, 
I13,i67,]j this is irrelevant, because the gluon loops are still 
IR sub-leading. 

It remains to select for the numerical solution at in- 
termediate momenta the Matsubara frequency for which 
the dressing functions should be evaluated in the function 
A. The quantity \/ql+pl + (go + PoY / {^t^T) is, apart 
from po = 90 = (the most relevant case), in general 
not an integer. As there is no obvious possibility, the 
largest integer smaller than this expression will be taken 
as the Matsubara frequency at which the expression is 
evaluated. 

Note that in the present approach the use of the Mat- 
subara formalism and of the vacuum perturbation theory 
in an asymptotically-free theory does not contradict the 
Narnhofcr-Thirring theorem [68|. In fact, solving the set 
of DSEs sclf-consistently is not equivalent to expanding 



23 



I Ghost dressing function | 




10"' 1 10 10" 

3d-transverse gluon dressing function | 



% 1.6 












N 1.4 






1.2 







1 






0.8 










T=10Aqcd.P„=0 




0.6 




.... T=0.25 Aqcd, P„=0 




0.4 




T=10Aqcd.P„=27iT 




0.2 




.... T=0.25 Aqcd, P„=2t 1 






/ 1 1 1 







10'' 1 10 10^ 

3d-longitudinal gluon dressing function 



10' 




FIG. 14: Numerical results from the Dyson-Schwinger equa- 
tions at two temperatures in units of Aqqd for the soft mode 
po = and for the first hard mode po = 27rT (4 and 5 modes 
have been calculated explicitly in the high and in the low 
temperature case, respectively). We show in the top, mid- 
dle and bottom panel the dressing functions of the ghost, 
of the 3d-transverse gluon and of the 3d-longitudinal gluon, 
respectively. The different ultraviolet behavior for different 
temperatures is due to the temperature-dependent choice of 
the renormalization conditions. 



around a free system and the obtained propagators thus 
describe thermal quasi-particles. 



2. Numerical solution 

A full numerical solution of the DSEs (for an explicit 
form of the equations see [33j ) and the present vertex con- 



struction can be performed using the method described 
in [g^. Let us note that, in order to improve the speed 
of the algorithm, it is useful to let the damping con- 
stants, introduced in Ref. ^69|, decrease with the itera- 
tion number, starting from a very large value. In prac- 
tice, however, only a (small) finite set of the Matsubara 
frequencies can be treated independently, due to limita- 
tions in computing power. For the Matsubara frequen- 
cies not determined explicitly we use the perturbative 
behavior reported in Eqs. (jD6HD7|) . Thus, at sufRciently 
large momenta, perturbation theory becomes dominant 
and the 0(4)-invariance is restored. We also approxi- 
mate the summation as an integration, starting from a 
given frequency (larger than the largest one determined 
independently). This integration extends to negative and 
positive infinity, respectively, and can be treated using a 
normal Gauss-Legendre integration. 



Results for two different temperatures are shown in 
Fig. [H for the S'C/(3) case, with ^ = 1 and C = 3 (i.e. 
K = —1/2) il3|]. These results are qualitatively similar 
to our findings using lattice calculations (see Section [1111) 
and to results reported in Refs. [3l|. A further inter- 
esting observation is the temperature dependence of the 
IR coefficients, e.g. the one of the ghost propagator (oq) 
increases with decreasing temperature. 



Note that the temperatures considered here are quite 
large. While there is no problem in going to arbitrar- 
ily large temperatures and obtaining the infinite temper- 
ature limit explicitly, going to smaller temperatures is 
limited by two technical problems. One is the sheer com- 
puting power and the necessity to include more Matsub- 
ara frequencies. The other one is a truncation-induced 
problem. Indeed, the mid-momentum behavior of the 
three-gluon vertex destabilizes the system, as the only 
possible solution would have a sign change in the gluon 
propagator, which is not allowed due to Eqs. In 
order to reduce this problem, the numerical calculations 
have been done using an interpolation between the per- 
turbative and the IR behavior, instead of considering the 
full dependence on the dressing functions. Such a change 
does not induce an error larger than the one already 
present due to the truncation scheme. On the other hand, 
this artifact can only be cured by an adequately chosen 
effective vertex-dressing with a stronger mid-momentum 
suppression. This is actually what is expected from re- 
cent lattice calculations 0, [S^, 113] . The same problem 
appears when considering ^ ^ 3. Thus, these results 
should be taken to be a proof-of-principle that this sys- 
tem indeed has solutions of the described type, rather 
than a quantitative investigation. 
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